Bose-glass phases of ultracold atoms due to cavity backaction 
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We determine the quantum ground-state properties of ultracold bosonic atoms interacting with 
the mode of a high- finesse resonator. The atoms are confined by an external optical lattice, whose 
period is incommensurate with the cavity mode wave length, and are driven by a transverse laser, 
which is resonant with the cavity mode. While for point like atoms photon scattering into the cavity 
is suppressed, for sufficiently strong lasers quantum fluctuations can support the build-up of an 
intracavity field, which in turn amplifies quantum fluctuations. The dynamics is described by a 
Bose-Hubbard model where the coefficients due to the cavity field depend on the atomic density 
at all lattice sites. Quantum Monte Carlo simulations and mean- field calculations show that for 
large parameter regions cavity backaction forces the atoms into clusters with a checkerboard density 
distribution. Here, the ground state lacks superfluidity and possesses finite compressibility, typical 
of a Bose-glass. This system constitutes a novel setting where quantum fluctuations give rise to 
effects usually associated with disorder. 
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Bragg diffraction is a manifestation of the wave- 
properties of light and a powerful probe of the micro- 
scopic structure of a medium: Bragg peaks are intrin- 
sically related to the existence of spatial order of the 
scatterers composing a medium and provide a criterion 
for the existence of long-range order |l|. Bragg diffrac- 
tion of light by atoms in optical lattices has been mea- 
sured for various geometries and settings, from gratings 
of laser-cooled atoms [sMHIl to ultracold bosons in the 
Mott-Insulator (MI) phase [(j. In most of these setups 
the backaction of light on the atomic medium, due to the 
mechanical effects of atom-photon interactions, is usu- 
ally negligible, while photon recoil can give rise to visible 
effects in the spectrum of the diffracted light Q. 

Recent work proposed to use high-finesse optical res- 
onators to enhance light scattering into one spatial di- 
rection, increasing the collection efficiency and thereby 
suppressing diffusion related to photon scattering [§| . For 
appropriate geometries, properties of the medium's quan- 
tum state can be revealed by measuring the light at the 
cavity output @, Q . These proposals assume that back- 
action of the cavity field on the atoms can be discarded. 
Such an assumption is, however, not valid in the regime 
considered in Refs. jol— Fl^U| : Here, the strong coupling 
between cavity and atoms can induce the formation of 
stable Bragg gratings in cold 0, [loj and ultracold atomic 
gases pj]-[l4[ that coherently scatter light from a trans- 
verse laser into the cavity mode. This phenomenon oc- 
curs when the intensity of the pump exceeds a certain 
threshold [l(| [HI, EH, [THj] . At ultralow temperatures the 
self-organized medium is a supersolid [13| , while for larger 
pump intensities incompressible phases are expected [16j . 

Let us now assume that the atoms are inside a high- 
finesse standing-wave cavity and form a periodic struc- 
ture, like the one sketched in Fig. [TJ whose period is 
incommensurate with the cavity mode wavelength. If 
the atoms are point like scatterers, or deep in a MI phase 




FIG. 1. (color online) Ultracold atoms are confined by an op- 
tical lattice of periodicity Ao/2. They are transversally driven 
by a laser at Rabi frequency Q and strongly coupled to the 
standing-wave mode of a high-finesse cavity with wavelength 
A. The wave lengths A and Ao are incommensurate: The 
light scattered into the cavity mode gives rise to an incom- 
mensurate potential, which mediates an effective long-range 
interaction between the atoms and can substantially modify 
the quantum ground state. 



of the external potential, there is no coherent scatter- 
ing from a transverse laser into the cavity mode [H, [I?} . 
Quantum fluctuations, however, will induce scattering 
into the cavity mode, thus the creation of a weak periodic 
potential which is incommensurate with the periodicity 
of the optical lattice confining the atoms. In this Letter 
we derive a Bose-Hubbard model for the system in Fig. 
[U and show that cavity backaction gives rise to an addi- 
tional term in the atomic energy which depends on the 
density at all lattice sites. Even for weak probe fields this 
results in the formation of patterns which maximize scat- 
tering into the cavity mode and can exhibit finite com- 
pressibility with no long-range coherence. This feature, 
typical of disordered systems, corresponds to a Bose-glass 
(BG) phase for sufficiently deep potentials [18-21] and 
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here emerges due to the nonlocal quantum potential of 
the cavity field. These predictions should be observable 
in existing experimental setups [l3|, [Hj . These findings 
extend studies where glassiness was predicted in multi- 
mode cavities [23|, |24[ . 

To analyze the quantum dynamics of the atom-cavity 
system we assume that the atomic motion is confined 
to the x — z plane where r = (x, z) is the atomic po- 
sition. It is furthermore assumed that the scattering 
processes are coherent, namely, the modulus of the de- 
tuning A a = ujl — ujq between the frequencies of pump, 
uol, and atomic transition, cjo, is much larger than (i) 
the transition linewidth, (ii) the strength of the atom- 
photon coupling, and (hi) the modulus of the detuning 
5 C = ojl — oj c between pump and cavity mode [l3|, fl6j . 
The laser and cavity mode have wave vector = kx 
and k cav = kz, respectively: When the polarizations are 
suitably chosen, only one stable electronic state is oc- 
cupied at all times. We denote by a and a) the oper- 
ators annihilating and creating a cavity photon and by 
^(r) the bosonic field operator for an atom at r, with 
|/0(r), ^(r')] = S^ 2 \r — r'). In the reference frame ro- 
tating at frequency ljl the Hamiltonian governing the 
coherent dynamics reads H = —h5 c a)a + H a + Hi, with 
Hi =a ^i = j dr0t(r)7^(r) , 0(r). Here, 7i a (r) describes the 
atomic dynamics in absence of the resonator: 

Ha(r) = -|^V 2 + Vb(r) + Q s n(r) + V 1 cos 2 (kx) , 
2m 

where Vb(r) = V t {cos 2 (koz) + f3 cos 2 (kox)} is an external 
periodic potential of wave number ko, depth V t along z 
and aspect ratio /3, n(r) = ^(r)^^) is the atomic den- 
sity operator, Q s is the strength of s-wave collisions, and 
V\ = Ml 2 j A a such that \V±\ is the depth of the dynam- 
ical Stark shift induced by the transverse standing-wave 
laser at Rabi frequency ft. The coupling with the cavity 
gives Hi(r) with jl6l l25j 



Hi(r) = HUq cos 2 (kz)a)a + HSo cos(kx) cos(kz) (a 



(1) 



where Uo = g^j A a is the dynamical Stark shift per cav- 
ity photon and go is the vacuum Rabi frequency at a 
cavity-field maximum. The frequency So = goft/A a is 
the amplitude of scattering a laser photon into the cav- 
ity mode by one atom. The corresponding term describes 
the coherent pump on the cavity field via scattering by 
the atoms and depends on the atomic positions within the 
field. It gives rise to a large intracavity-photon number 
n cav when the atoms form a Bragg grating with period- 
icity 27r/fc. This would correspond to choosing k = ko 
[HI : n cav would then depend on the balance between the 
superradiant scattering strength, SoN, and the rate of 
cavity loss n. On the contrary, in this paper we take k 
and ko to be incommensurate and analyse the effect of 
the strong coupling with the cavity field when the atoms 
are tightly confined by the potential Vb(r) and weakly 
pumped by the transverse laser. 



We now sketch the derivation of an effective Bose- 
Hubbard Hamiltonian describing the dynamics of the 
atoms in the cavity potential. We first assume that 
in the time scale At the atomic motion does not sig- 
nificantly evolve while the cavity field has reached a 
local steady state, J^ At d(r)dr/ At w a st . Here, 

J^ +A£ d st (r)dr = where the time evolution is gov- 
erned by the Heisenberg-Langevin equation of motion 
a = [d,H]/ih — kcl + V^&in^), with d[ n (t) the input 
noise such that [d[ n (t), d} n (t f )] = S(t — t'). The condi- 
tions for a time scale separation are set by the inequal- 
ity | S c + At ^> 1 and by the assumption that coupling 
strength between atoms and fields are much smaller than 
1/At [27]. The stationary field reads 



Oat 



SoZ 



iv2^ai 1 



(s c - u y) + in (s c - u y) + 



(2) 



with a- in the input noise averaged over At. Here, y = 
J d 2 r cos 2 (kz)h(r) and Z = f d 2 r cos(kz) cos(kx)h(r) 
describe the shift of the cavity resonance and coher- 
ent scattering amplitude, respectively, due to the atomic 
density distributions. The quantum noise term can be 
neglected when the mean intracavity-photon number is 
larger than its fluctuations, that corresponds to taking 
\So(Z}\ ^> n. In this limit, the field at the cavity out- 
put, a ou t = V^^st — ain, allows one to non-destructively 
monitoring the state of the atoms [I?}- Using Eq. 
in place of the field a in Eq. (pQ) leads to the effec- 
tive atomic Hamiltonian, from which we derive a Bose- 
Hubbard model assuming that the atoms are tightly 
trapped in the lowest band of Vb(r). We use the Wannier 
decomposition ^(r) = Wi(r) hi, with ^(r) the lowest- 
band Wannier function centered at the minimum v\ of the 
classical potential Vb (r) and hi the bosonic operator anni- 
hilating a particle at n, such that hi = b\bi is the onsite 
atomic density. Using the thermodynamic limit where 
the cavity parameters scale with the number of lattice 
sites K according to So = so/VK and Uo = uo/K [liSj . 
the resulting Bose-Hubbard Hamiltonian reads 



U , 



Hbh = ^2 \~ ^2 iiffl-bj + H.c.) + ^hi(hi - 1) - forii 

(3) 

with U = Q s J d 2 r [wi(r)] A the onsite interaction 
strength, fa the site-dependent chemical potential, 
U the tunneling coefficient, and the sum over 

nearest neighbors. Here, fa = /i^ ^ + Sfii, ii = 
+ Sii, where ^ = -E - V t X and t^ = 
—Ei — V t Xi are constant over the lattice, with 
E s=0 ,i = -h 2 /(2m) J d 2 r^(r)V 2 ^ i+s (r) and X s=0 ,i = 
J d 2 r {cos 2 (koz) + f3 cos 2 (kox)}wi(r)wi+ s (r). The terms 
5 ^ and 5U are due to the pump and cavity incommen- 
surate potentials and vanish when the pump laser is off, 
Q = 0. In this limit the model reduces to the regu- 
lar Bose-Hubbard model exhibiting a Ml-superfluid (SF) 
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FIG. 2. (color online) Results of QMC simulations for a ID lattice (f3 ^> 1) with 100 sites and periodic boundary conditions, (a) 
Phase diagram in the ji — t plane for sq = 0.004ft. The gray regions indicate incompressible phases at density n = 1,2, the blue 
region the gapless phases with vanishing SF density. The dotted lines indicate the shape of the incompressible phases in absence 
of the pump laser (so = 0) [3(J] . (b) Linear density h and (c) (i>) as a function of \i for t — 0.053£7 and so/k = 0.003 (triangles), 
so/k = 0.004 (circles), and so = (squares): The number of photons is different from zero for the parameters corresponding to 
the blue regions in (a). The inset of (b) displays the Fourier transform of the pseudo current-current correlation function J[oS) 
[3(3] for the parameters indicated by the arrows in the curve of (b). (d) Local density distribution (hi) (empty points joined 
by the blue curve) and local density fluctuations (n?) — {hi) 2 (filled points joined by the red curve) as a function of the site 
for fi — and so = 0.004ft. The values of so are found from so = VKQgo/ A a and are consistent with the parameters of Ref. 
O (po/2tt = 14.1 MHz, A a /27r = 58 GHz, ft/27r = 1.3 MHz). The other parameters are 6 C = -5ft, k/k = 785/830 as in 0, 
while Q s has been taken from Ref. [33l ]. 



transition as a function of the parameters fi and t [18| . 
When Q ^ 0, instead, one has 

Sfc = - Vl jW - hj^-^ (2<5 eff Zf + u l>lf ) , (4) 

"eff + K 

while oU is neglig ible, hence U = [28]. In Eq. 
(jH we introduced = J d 2 r cos 2 (kx)wi(r) 2 , which 

scales the strength of the classical transverse potential 
due to the laser (it is constant for the sites with same x 

value). The coefficients for the cavity field read Z$ = 

/ d 2 r cos(kx) cos(kz)wi(r) 2 , = J d 2 r cos 2 (kz)wi(r) 2 , 

while operator <5 e ff = S c — uq hi/ K accounts for 

the cavity- frequency shift due to the atoms [25|, [29| . Re- 
markably, the cavity effects are scaled by the operator 

K 

$ = '£4 i) ni/K, (5) 

i=l 

which originates from the long-range interaction medi- 
ated by the cavity field, and is related to the mean in- 
tracavity photon number since n cav oc (4 )2 ). Its mean 
value vanishes when the atomic gas forms a MI state: 
(4 > )mi oc Zq^ = 0, since there is no coherent scat- 
tering into the cavity mode. Also deep in the SF phase 
($)sf —> 0- Close to the MI-SF phase transition, how- 
ever, fluctuations in the atomic density lead to finite val- 
ues of (4>), hence to a finite intracavity photon number. 
The dependence of the chemical potential on the operator 
<l> is a peculiar property of our model, that distinguishes 
it from the case of a bichromatic optical lattice [20|, Hl| , 
where the strength of the incommensurate potential is 
independent of the phase of the ultracold atomic gas. 



We first analyse the ground state of the Bose-Hubbard 
Hamiltonian, Eq. (|3j), when the system can be reduced 
to a one-dimensional (ID) lattice along the cavity axis 
(f3 ^> 1). In ID the effect of the transverse potential 
(first term in Eq. @) is a constant shift which can be 
reabsorbed in the chemical potential. Figure [2fa) dis- 
plays the phase diagram in the n — t plane evaluated 
by means of a quantum Monte Carlo (QMC) calculation 
|30j . The compressibility is determined using \ = 
with n = ^2i(ni)/K the mean density, while the SF den- 
sity is obtained by extrapolating the Fourier transform 
of the pseudo current-current correlation function J(uo) 
[30| to zero frequency (see inset of panel (b)). The grey 
regions indicate the MI states at densities n = 1,2, the 
blue regions a compressible phase where the SF density 
vanishes, while outside the phase is SF. The effect of cav- 
ity backaction is evident at low tunneling, where (<£) > 0: 
Here the size of the MI regions is reduced. At larger tun- 
neling a direct MI-SF transition occurs and the MI-SF 
phase boundary merges with the one found for so = 0: 
In fact, for larger quantum fluctuations (4>) — >• 0. This 
feature is strikingly different from the situation in which 
the incommensurate potential is classical [34],[35|: There, 
the MI lobes shrink at all values of t with respect to the 
pure case. The atomic density and corresponding value 
of (4>) are displayed in Fig. EJb) and (c) as a function 
of fi for different values of the transverse laser inten- 
sity (thus so): the incommensurate potential builds in 
the BG phase of the diagram, where the number of in- 
tracavity photons does not vanish. Panel (d) displays 
the local density distribution (hi) and the local density 
fluctuations in the BG phase: the density oscillates in a 
quasi-periodic way over clusters in which the atoms scat- 




FIG. 3. (color online) Results of a mean-field calculation for a 2D lattice with 70 x 70 sites and periodic boundary conditions, 
(a) Order parameter in the \i — t plane when so = 0.15ft (S c = —5ft). Inset: phase diagram when cavity backaction is 
negligible (S c — —300ft) 39]. The dotted (solid) curve delimitates the regions with a small, arbitrarily chosen threshold for 
the order parameter (density fluctuations) [40]. The dashed curve shows the region where the number of photons is 2 orders 
of magnitude larger than outside, (b) Local density distribution (hi) at ji = 0.156?7 and t = 0.01U (point in (a) at n = 0.64). 
Inset: corresponding density distribution in the ji — t plane when cavity backaction is negligible (here n = 0.625). (c) Photon 
emission rate at the cavity output, n ou t = (aJ ut a ou t), as a function of t for n = 0.64. The rate is in units of n^ t = ftnc^ ax \ 
where ricSf 3 ^ = SqK/(5^ + ft 2 ) is the maximum number of intracavity photons when all atoms scatter in phase into the cavity 
mode. Inset: corresponding order parameter. Note that in absence of cavity backaction (5 C = —300ft) the rate is less than six 
orders of magnitude smaller. The value so — 0.15ft is consistent with the parameters of Ref. [l3|] when taking 300 x 300 lattice 
sites. The onsite repulsion Q s is taken from Ref. [4ll ]. 



ter in phase into the cavity mode. The fluctuations are 
larger at the points where Zq \ which oscillates at the 
cavity mode wave length, becomes out of phase with the 
trapping potential. In this way the density distribution 
maximizes scattering into the cavity mode. 

The two-dimensional ( 2D) case is studied by means of 
a mean-field calculation [l8|, [36| for a pumping strength 
so such that max{|(5/ii)|} < U. Here, the transverse 
laser determines the transverse density distribution even 
when the cavity field is zero. In this case the MI lobes at 
t = shrink due to the classical incommensurate field. 
Due to the classical potential, in general no direct MI- 
SF transition is expected [37j . Figure [3fa) displays the 
order parameter, 52i(bj)/K, in the /i — t plane and for 
density n < 1 [38|. The solid curve indicates where the 
gap in the spectrum is different from zero, correspond- 
ing to vanishing density fluctuations A^> = (n 2 — n 2 ) 1 / 2 , 
for n 2 — ^i(n 2 ) / K [42|, |43|- The dashed line separates 
the parameter region where the number of intracavity 
photons is at least 2 orders of magnitude larger than 
outside. For comparison the diagram for n cav = is re- 
ported in the inset. A larger region with vanishing order 
parameter appears when the atoms are strongly coupled 
to the cavity. The analysis of the corresponding density 
distribution shows a distinctively behaviour. Panel (b) 
displays the local density for 5 C = —5ft for the parame- 
ters indicated by the point in (a): The coupling with the 
resonator induces the formation of clusters with checker- 
board density distribution, where the atoms scatter in 
phase into the cavity mode. At the border of the clusters 
the density fluctuates, so to allow the fields scattered by 
each cluster to interfere constructively. The inset shows 



that the clustering disappears when cavity backaction is 
negligible. In this latter case the field at the cavity out- 
put is zero, while when the effect of cavity backaction 
is relevant the corresponding field at the cavity output 
(panel (c)) has a finite intensity for a finite range of val- 
ues of t: When the tunneling rate is instead sufficiently 
large, there is no clustering and the atomic density be- 
comes uniform along the cavity axis. Thus, the compe- 
tition between the mechanical effects of the cavity field 
and quantum fluctuations leads to the creation of these 
clusters, which exhibit properties ranging from a BG to 
a strongly-correlated SF, and then disappear when tun- 
neling becomes large. 

These results shed new light onto the effect of cav- 
ity backaction on an atomic ordered structure. At ul- 
tralow temperatures, when the ordered medium does not 
support photon scattering into the cavity mode, cavity 
backaction forces the system into locally ordered clusters 
which are phase locked with one another, thereby increas- 
ing the intracavity photon number. This is an example 
of selforganization of the quantum gas that is triggered 
by the quantum fluctuations of the atomic motion. 
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